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We model the dynamics of social structure by a simple interacting particle system. The social 
standing of an individual agent is represented by an integer-valued fitness that changes via two 
offsetting processes. When two agents interact one advances: the fitter with probability p and the 
less fit with probability 1 — p. The fitness of an agent may also decline with rate r. From a scaling 
analysis of the underlying master equations for the fitness distribution of the population, we find 
four distinct social structures as a function of the governing parameters p and r. These include: (i) 
a static lower-class society where all agents have finite fitness; (ii) an upwardly-mobile middle-class 
society; (iii) a hierarchical society where a finite fraction of the population belongs to a middle class 
and a complementary fraction to the lower class; (iv) an egalitarian society where all agents are 
upwardly mobile and have nearly the same fitness. We determine the basic features of the fitness 
distributions in these four phases. 
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I. INTRODUCTION 

The emergence of class structure in society is a ubiq- 
uitous phenomenon in the biological and the social sci- 
ences 0, 0, S 01 • Social hierarchies have been widely ob- 
served in animal populations including insects [f| , mam- 
mals 0, 0J3 , and primates , as well as human com- 
munities |lfT| . 

The possibility of quantitative modeling of social phe- 
nomena using concepts and techniques borrowed from 
the physical sciences is rapidly gaining appreciation. Ex- 
amples of such modeling include the wealth distribution 
op inion dynamics [l3L IfH Il5| , and rumor propa- 
gation |16lll7| . Such approaches typically draw analogies 
between individual agents in the social system and par- 
ticles in a corresponding physical system and then iden- 
tifying macroscopically observed phenomena with micro- 
scopic agent-agent interactions 0, 0, |2(| • 

In this spirit, we seek to tie the emergence of social 
structures to specific interactions between agents within a 
general version of the recently-introduced advancement- 
decline process 0, In our model, the social 
standing of each agent is characterized by a single num- 
ber, its fitness. Agents increase their fitness by interact- 
ing with other agents and also, their fitness may decline 
spontaneously. This simple model has only two parame- 
ters: the probability that the fitter agent advances in an 
interaction and the rate of decline. 

We find that a rich variety of familiar social structures 
emerges as a result of the competition between advance- 
ment and decline. When decline dominates, the society 
is static and the fitness distribution approaches a steady 



state. When the decline rate is comparable to the ad- 
vancement rate, the society is dynamic and the charac- 
teristic fitness of the population increases linearly with 
time. In this case, there are several possibilities. When 
the less fit agent benefits from social interactions, an egal- 
itarian society arises in which all agents advance at the 
same rate. Consequently, inequalities among agents are 
small. On the other hand, when the fitter agent tends to 
benefit in competitions, agents advance at different rates 
and social inequalities increase with time. Depending on 
the relative influence of advancement and decline, either 
the entire population or only a fraction of it may be up- 
wardly mobile. In the latter case, the society consists of 
a static lower class and an upwardly-mobile middle class. 

In Section ^ we introduce the general advancement- 
decline process and the governing master equations. The 
overall class structure and the statistics of the mobile 
middle class are obtained using scaling analysis in Sec- 
tion [ffi] The basic features of the egalitarian society are 
investigated in section llVl where the cumulative fitness 
distribution may be largely determined by linear trav- 
eling wave analysis. In Section the statistics of the 
lower class, where the fitness distribution is steady, are 
determined. We conclude in section IvH 
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II. THE ADVANCEMENT-DECLINE MODEL 

We model a scenario in which the social status of an 
agent benefits from increased social interactions, while 
solitude or isolation have the opposite effect. Indeed, 
highly connected individuals often have better access to 
information, resources, and power, that are often gained 
as a result of social interactions. Thus, in our model there 
are two competing evolutionary processes that influence 
the fitness of agents: (i) advancement via social inter- 
actions, and (ii) decline due to the lack of interactions 
(Fig. ^| . For simplicity, social standing is represented by 
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a single parameter, the integer-valued fitness k > 0. 

(i) Advancement. Agents interact in pairs, and as a re- 
sult of the interaction, only one advances. There are two 
possibilities: either the fitter agent advances or the less 
fit advances. We allow the fitter agent to advance with 
probability p and the less fit agent to advance with prob- 
ability 1 — p. Thus, when two agents with fitness k and 
fitness j interact (with k > j), the outcome is 

— * (k + X,j) with probability p, 
(k,j) — > (k,j + f) with probability 1 — p. 

For p — 1 the fitter agent always advances |2^| , while for 
p = the less fit agent always advances. The interac- 
tion rules are defined so that one randomly-chosen agent 
advances when two equally-fit agents interact. Without 
loss of generality, the interaction rate is set to 1/2. Also, 
we consider the thermodynamic limit where the number 
of agents is infinite. 

(ii) Decline. In the decline step, the fitness of an indi- 
vidual decreases according to 

fc->fc-l 

with rate r. This process reflects a natural tendency for 
fitness to decrease in the absence of social activity. We 
impose the lower limit for fitness to be k = 0; once an 
individual reaches zero fitness, there is no further decline. 

Our goal is to understand how the fitness distribution 
of a population evolves as a function of the two model pa- 
rameters, the advancement probability p and the decline 
rate r. Let fk(t) be the fraction of agents with fitness 
k at time t. In the mean-field limit, where any pair of 
agents is equally likely to interact, the fitness distribution 
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obeys the master equation 

-77 = K/fc+i - fk) +P(fk-iFk-i - fkF k ) 
dt i (1) 

+ (1 - p)(fk-xGk-x - f k G k ) + (fk-i - fk) ■ 

Here F k = J2jZo fj and G k = J^JLk+i fo are tne respec- 
tive cumulative distributions of agents with fitness less 
than k and fitness greater than k. The boundary condi- 
tion is f-i(t) = 0. The first pair of terms accounts for 
decline, the second pair of terms describes interactions 
where the stronger agent advances, and the third pair of 
terms accounts for interactions where the weaker agent 
advances. The last pair of terms describes interactions 
between two equal agents and it reflects that when two 
such agents interact, only one of them advances. The 
prefactor 1/2 arises because there are half as many ways 
to chose equal agents as there are for different agents. 
We consider the initial condition where all agents have 
the minimal fitness /fc(0) = 5k t o- 

It proves useful to rewrite the evolution equation in 
a closed form that involves only the cumulative distri- 
bution. Summing the rate equations Q and using the 
relations fk = -Ffe+i — F k and Gk = 1 — F k+ \, the cumu- 
lative distribution F k obeys 

dFi 

-rf = r(F k+ i - F k ) + pF k _ 1 {F k - 1 - F k ) 
dt (2) 

+ (1 - p)(l - F k ){F k _i - F k ) - -{F k - iVx) 2 . 

The boundary conditions are Fq — 0, Foo — 1, and the 
initial condition is F k (Q) = 1 for k > 1. There is a one-to- 
one correspondence between the four terms in equations 
JTJ and J2J • The master equation for the cumulative dis- 
tribution can be simplified by consolidating the advance- 
ment terms 

^=r(F k+1 -F k ) + (l-p)(F k _ 1 -F k ) 



k-2 k-1 k k+1 k+2 + (p - 1/2) {F k _ x - Fl) . 
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The mean fitness (k) = ^\ kf k evolves according to 




FIG. 1: The elemental processes of the advancement-decline 
model. Top: the decline process. Bottom: the advancement 
process, with either the fitter agent advancing (right) or the 
less fit agent advancing (left). 



This result can be derived directly by summing the mas- 
ter equations Q or even simpler, from the definition 
of the advancement-decline process. The first term ac- 
counts for advancement, where interactions occur with 
rate 1/2 such that each interaction advances only one 
agent. The second term stems from decline and reflects 
the fact that all agents except for the least-fit ones decline 
with rate r. 

We now discuss the basic social structures that emerge 
from the solution to the master equation. 



III. EMERGENCE OF SOCIAL STRUCTURES 
A. Scaling solution 

We determine the class structure of the population via 
a simple scaling analysis of the master equation. Let us 
take the continuum limit of the master equation by re- 
placing differences with derivatives, -Ffc+i — Fj. — > dF/dk. 
To first order in this "spatial" derivative, we obtain the 
nonlinear partial differential equation 

OF dF 

- = lp + r -l-(2 P - m -. (5) 

When the spatial derivative and the temporal derivative 
balance, the typical fitness increases linearly with time, 
k ~ t. Therefore, we make the scaling ansatz 



F k (t) ~ $ 



(6) 



The boundary conditions are $(0) = and $(00) = 1. 

Substituting this scaling form in Eq. the partial- 
differential equation reduces to the ordinary differential 
equation 

[(p + r-l + 3!)-(2p-l)*(»)] — =0, (7) 

where the prime denotes differentiation with respect to 
the scaling variable x = k/t. The solution is cither 
d®/dx = 0, i.e., 



$(2;) = constant, 



or the linear function 



2p - 1 2p - 1 



(8) 



(9) 



Using these two solutions and invoking (i) the bound- 
ary conditions ^(O) = and $(00) = 1, (ii) the bounds 
< <f>(x) < 1, (iii) monotonicity of the cumulative dis- 
tribution, d&(x)/dx > 0, and (iv) the assumption that 
the scaling function changes continuously with p and r, 
we can then deduce the four possible social structures of 
the population. 

1. Middle-Class Society: We first examine the conditions 
for the linear scaling function to apply. First, the cu- 
mulative scaling function @ must be a monotonically in- 
creasing function. Therefore, the linear solution holds 
only when its slope is positive, that is, when p > 1/2. 
Second, the scaling function is bounded, < $(2;) < 1; 
this condition implies the lower and upper bounds 



= 1 - (p + r) 



and 



x + =p- 



(10) 



on the scaled fitness. The obvious constraints X- > 
and x + > lead to the conditions p + r < 1 and p > r. 
By imposing continuity as well as the limiting behav- 
iors <&(&) = and <j>(x) = 1 outside the linear region, 
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FIG. 2: Phase diagram of the advancement-decline model. 
The small graphs in each region are sketches of the scaled 
cumulative fitness distribution. 
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FIG. 3: Middle-Class society. The scaled cumulative fitness 
distribution $(2;) versus x = k/t at different times for p = 5/8 
and r = 1/4. 



the scaled cumulative distribution is the piecewise linear 
function (Fig. EJ: 



$m{x) 





p + r — 1 x 
2p-l + 2p- 1 

1 



< x < x_ 

X- < x < x + (11) 

X + < X. 



This behavior describes a middle class society where all 
agents are upwardly mobile, as their fitness improves lin- 
early with time. In this case, social inequalities also in- 
crease indefinitely with time: the agents at the bottom 
of the middle class have fitness fc_ = [1 — (p + r)]t and 
the richest agents have fitness k + — (jp—r)t. The middle- 
class society lies within the triangular region defined by 
the lines r + p = 1, r = 0, and p = 1/2, shown in Fig. [21 

2. Hierarchical Society: Along the line r + p = 1, the 
fitness of the poorest agents vanishes. Moreover, the lin- 
ear scaling solution @ has a finite positive value at zero 
fitness for a range of parameter values p and r. These 
two observations suggest the existence of another type 
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FIG. 4: Hierarchical society. The scaled cumulative fitness 
distribution $(3;) versus x — k/t at different times for p = 3/4 
and r = 1/2. 



of piecewise linear solution with $(0) > 0. The bounds 
< $(0) < 1 impose the conditions p + r > 1 and r > p. 
In this region, the scaling function has two distinct com- 
ponents (Fig.0J 

!P + r — 1 x 
1 < x < x+ , x 
2p-l 2p-l + (12) 
1 x+ < x. 

Thus, we find a hierarchical society (Fig.© that includes 
both an upwardly-mobile middle class and a static lower 
class. The lower class consists of a finite fraction 



p + r — 1 
2p- 1 



(13) 



of agents with zero fitness (in scaled units). In sect ion Ivl 
we examine the lower class more closely and show that 
its fitness distribution is time-independent and extends 
only over a finite range. 

3. Lower- Class Society: When the fraction L of agents 
with zero fitness reaches 1, the entire population is poor. 
Forp > 1/2, the condition L = 1 occurs on the boundary 
p = r. At this point the fitness distribution becomes a 
step function, 



$ h (x) = e(x) 



(14) 



with 6(x) = for x < and G(x) = 1 for x > 0. We 
therefore conclude that there is a region of the phase di- 
agram where the scaled fitness of the entire population is 
zero. For any initial state, the fitness distribution quickly 
approaches the step-function in a lower-class society. 

4- Egalitarian Society: There is another region of the 
phase diagram where the fitness distribution also be- 
comes a step function. When p = 1/2 and r < 1/2, 
then both 2+ and x_ are equal to 1/2 — r. Therefore 
3>(x) = <d(x ~ v) with v = 1/2 — r. Since the scaling 
function must change in a continuous fashion, we con- 
clude that for p < 1/2, the scaling function is again a 



step function but with a jump at non-zero fitness. That 
is 



$e(x) = 0(x - v). 



(15) 



In this egalitarian society, all agents have the same scaled 
fitness x = v or alternatively the fitness k w vt. The ve- 
locity v follows easily from the average fitness (©. Since 
all agents advance at constant rate, then the term — r/o 
is negligible and therefore, the propagation velocity is 



(16) 



In section lTVl we show that in this society, the fitness dif- 
ferences between agents are small and do not grow with 
time. This is the sense in which the society is egalitarian. 
When p < 1/2, the weaker agent preferentially benefits 
in an interaction, so that the rich effectively supports the 
poor. We also note that the lower class and the egalitar- 
ian society share one common feature: they do not have 
a middle class. The boundary between these two phases, 
determined by the condition v = 0, is the line r = 1/2 
(Fig.©. 

Our numerical integration of the evolution equations 
confirms the overall picture of four different social struc- 
tures (Fig. © : a middle class society (Fig. © , a hierarchi- 
cal society (Fig.©, a lower-class society as in ifHj) . and 
an egalitarian society as in l|15l) . The numerical data was 
obtained by integrating Fk for < k < 10000 using a 
fourth-order Adams-Bashforth method with a spec- 
ified accuracy of 10~ 10 in the distribution Fk- 



B. Refinements to the Scaling Solutions 

Our numerical results for the cumulative distribution 
Fk, when plotted versus the scaling variable x = k/t, 
smoothly approaches the appropriate expressions for the 
piecewise linear scaling function $(x) derived in the pre- 
vious subsection (Figs. |3] & QJ. As time increases, the 
fitness distribution narrows. The simulations also show 
that the approach to the scaling solution is slowest in the 
vicinity of the extremes of the middle class x — X- and 
x = x + (with x_ = for the hierarchical phase). 

The correction to scaling near these extrema can be 
determined by keeping derivatives up to second order in 
the continuum limit of the master equation. This approx- 
imation gives the nonlinear diffusion equation l2fj| | 



dF 
~dt 



P) 



d 2 F 
~Wk 



(17) 



(1 - 2p)F 



dF _ l d 2 F 
~dk ~ YWk 



+ (p-l/2) 



dF 
~dk 



The linear terms are separately displayed in the first line 
and the nonlinear terms in the second. 

Let us first consider the poorest agents, i.e, the be- 
havior close to x = X-. Since the cumulative fitness 
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distribution is small near this point, the nonlinear terms 
can be neglected and the governing equation l(T7|l reduces 
to the standard convection-diffusion equation 



dF 
~dt 



V- 



dF 
~dk 



d 2 F 



(18) 



with propagation velocity u_ = 1 — p — r and 

diffusion coefficient £)_ = (1 — p + r)/2. Indeed, since 
the fitness distribution is obtained from the cumulative 
distribution by differentiation, / = dF/dk, the fitness 
distribution satisfies the same equation JSJ as the cumu- 
lative distribution. 

For the middle-class society, we therefore conclude that 
the bottom of the middle class has a Gaussian tail, with 
the center of the Gaussian located at fc_ = v-t and with 
width J D-t. The same analysis can be carried out for 
the hierarchical society, where the quantity F — L now 
satisfies the diffusion equation with zero velocity u_ = 
and diffusivity D_ = r. Conversely, the distribution for 
the top end of the middle class can be obtained by ana- 
lyzing I — F. It is immediate to show that this quantity 
again obeys Eq. I|18|) with velocity v + = x+ = p — r 
and diffusivity D+ = (r + p)/2. We conclude that the 
extremes of the middle class are characterized by Gaus- 
sian tails whose extents grow diffusively with time. In 
terms of the scaling variable x, the deviation from the 
scaling function <&(x) is appreciable only within a region 
of whose width is shrinking as t^ 1 ! 2 . 

For the special case p — 1/2, the nonlinear terms van- 
ish and the fitness distribution is described exactly by the 
linear convection-diffusion equation l|I8l) with drift veloc- 
ity v = 1/2 — r and diffusion coefficient D = (r +p)/2 
(the nonlinear term is negligible). Thus there is a drift 
toward smaller fitness for r > 1/2 and the fitness dis- 
tribution approaches a steady-state profile that decays 
exponentially with fitness. In the opposite case of v > 0, 
the fitness distribution is simply a Gaussian that drifts 
to larger fitness with velocity \ — t and whose width is 
proportional to \fT5i. In the case of p = 1/2, the rela- 
tive position of an agent in the society is irrelevant and 
advancement reduces to a pure random walk |27| . 



IV. EGALITARIAN SOCIETY 

In the egalitarian phase, the step function form of the 
scaling solution, Eq. l(T5|) . suggests that the fitness dis- 
tribution has the traveling wave form 



F k (t)^U(k-vt) 



(19) 



with the propagation velocity l|16[) . This is confirmed 
by numerical integration of the master equation J3J), as 
shown in Fig. (J5J. 

To determine the shape of the wave U (z) analytically, 
we substitute the waveform ()I9|) into the master equation 
© to give the nonlinear difference-differential equation 
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FIG. 5: The scaled cumulative fitness distribution U (z) versus 
z = k — vt, where v = 1/2 — r is the speed of wavefront at 
different times for p = 1/4 and r = 1/4 (egalitarian society) . 



for U{z) 

-vU'(z) = r[U(z+l)-U(z)] + (l-p)[U(z-l)-U(z)} 
+ (p-l/2)[U 2 (z-l)-U 2 (z)]. (20) 

The boundary conditions are U(—oo) = and U(oa) = 1. 



A. Waveforms in the tail regions 

We apply standard linear analysis in the tail regions to 
deduce the leading and trailing shapes of the waveform. 
When z — > — oo, then U<1 and therefore U 2 -C U. To 
first order in U, Eq. I|20|l becomes 

vU' + r[U(z + l) - U(z)} + (l-p)[U(z-l) - U(z)} = 0. 

The behavior in this case is determined by the balance 
between decline and advancement events where the less 
fit agent advances. The solution to this linearized equa- 
tion is the exponential decay 



U(z) 



-oo. 



(21) 



Substituting this form and (|16|) into the linearized equa- 
tion, the decay constant a is the root of the following 
equation 

1/2 - r = a' 1 [(1 - p)(l - e- a ) - r(e a - 1)] . (22) 

Similarly, in the limit z — > oo we linearize the wave 
equation (|2U|I for the small quantity R = 1 — U to obtain 

vR' = r[R(z) - R(z+1)\ + p[R{z) - R(z-1)\. 

In this case the behavior at large fitness is governed 
by the balance between decline and advancement events 
where the fitter agent advances [27j ■ The solution to the 
above differential equation is again the exponential decay 



R{z) 



-0Z 



oo, 



(23) 
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FIG. 6: Tails of the fitness distribution versus z = k — vt 
for p = r = 1/4. The theoretical predictions (I2H - 124H with 
a = 0.535572 and (3 = 0.930821 are indicated by straight 
lines. 



FIG. 7: The super-exponential tail: the quantity 2 z R(z) 
versus 2 Z , obtained by numerical integration of 112"^ to time 
t = 100 with r = 0. 



with the decay constant [3 satisfying 

1/2 - r = J3' 1 [p^ 13 - 1) - r (1 - e"?)] . (24) 

We conclude that the likelihood of having agents that 
are much richer or much poorer than the average fitness 
k — vt in the egalitarian society is exponentially small, 
as illustrated in Fig. [fJl The society therefore consists of 
agents whose fitnesses are all roughly the same, k ~ vt. 
As one might naturally anticipate, social inequalities are 
small under the dynamics in which the rich preferentially 
gives to the poor. 

B. Less fit advances (j> = 0) 

For the case where the less fit agent always advances, 
the fitness distribution has a special form. In this case, 
the complementary cumulative distribution obeys 

^ = r(G k+1 -G k ) ^-(GU-Gl) (25) 

with the initial condition G k (0) = 5 k) —i and the bound- 
ary condition G_i(i) = 1. 

We expect that the fitness distribution will continue to 
have the form of a propagating wave. Substituting the 
traveling wave form G k {t) — > R(k — vt) into the master 
equation (25} gives 

-vR'(z) = r [R{z + 1) - R(z)} + )- [R 2 {z - 1) - R 2 (z)] . 

An exponential solution does not give asymptotic balance 
of terms as z — > oo, and we therefore attempt a solution 
of the form R(z) ~ ip{z)e~^ z \ Substituting this form 
into the above equation and keeping only the dominant 
term ^R 2 (z — 1) on the right-hand-side gives 

viP{z)(f>'{z)e-^ « ^ 2 {z - l)e- 2 ^- x ). (26) 



For the positive terms on the left and the right hand side 
to balance, the dominant exponential terms must first 
balance, yielding the recursion equation <j)(z) = 2c/)(z— 1). 
The solution is the exponential <p(z) = C2 Z . Balancing 
the prefactors, vip(z)<p'(z) = \i> 2 {z — 1) yields tp(z) = 
81n2Ci>2 2 . As a result, the decay in the tail region is 
super-exponential j2^ 

R(z) - 81n2C*w 2 Z exp (-C*2 Z ) , (27) 

as z — > oo. The constants C and v should be determined 
numerically. Hence, the front of the traveling wave is 
extremely sharp. This tail characterizes statistics of the 
rich, so when the rich never benefits from interactions 
with the poor, rich agents are ultra-rare (Fig. EJ. Even 
though the leading tail extends to only a handful of sites, 
it is still possible to verify the super-exponential decay 
(|27fl. In contrast, the z — > — oo tail that characterizes the 
poor is not altered; it has the same exponential tail as in 
Eq. (23}. 



V. THE LOWER CLASS 

We now determine the fitness distribution of poorest 
agents, a class that exists in both the hierarchical and 
lower-class societies. As we shall now show, the fitness 
distribution of the lower class in the limit of small but 
non-zero fitness approaches a steady state. For the case 
of the hierarchical society, we write Fk = L{l—g k ), where 
L is the lower-class fraction, with the deviation gk van- 
ishing for large k. Substituting this form into the master 
equation (3J and setting the time derivative to zero, gives 

r 9k ~ 9h+1 =l-p+(p- 1/2)L (2-g k - g k _ x ) . (28) 
9k-i~ 9k 

Consider first the lower-class society, for which the en- 
tire population is poor, L = 1. Using this fact, and 
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neglecting terms of the order g 2 , we find 



9k - gfc+l 
9k-i - 9k 



P 

T 



(29) 



The solution to this equation is simply the exponential 
form g k ~ a k with a = p/r. Since f k — Fk+i — F k — 
9k - 9k+i, then 



(30) 



In the lower-class society, the fitness is confined to a very 
small range. Notice also that this exponential decay co- 
incides with the traveling wave solution l|23|l with v set 
equal to zero, as the decay function is now = p/r. 

Finally, we consider the hierarchical society. Using 113|) 
and following the same steps that led to Eq. (|29|l . we 
obtain 



9 k - 9k+i 
9k-i - 9k 



1 



19k 



(31) 



with 7 = (r + p — l)/r. To determine gk, we expand 
the differences to second order and assume that g" <C 
g' to give, after straightforward steps, g" + jgg' = 0. 
The asymptotic solution to this equation is g ~ 2/(7/1). 
Finally, using f k = F k+1 - F k and F k = L(l-g k ), we 
determine the fitness distribution from fk — —Lg' to be 



fk 



2r 



2p-l 



(32) 



Thus, for the hierarchical society, the fitness distribution 
has a power-law large-fitness tail in the lower class region 
(see also for more details). 

As discussed in section llTTl there is a diffusive boundary 
layer that separates the steady-state fitness distribution 
in the lower class and time dependent fitness distribution 
in the middle class. From Eqs. © and i|12|) . the fitness 
distribution in the middle class is fk — [(2p — l)t] 
Equating this expression with Eq. ll-'il't gives a crossover 
scale 

K ^ V2ri. (33) 

Thus, the steady-state region extends over a fitness range 
that grows as t 1 ^ 2 . We also note that this crossover scale 
agrees with the diffusivity £>_ = r, obtained in section 
IIIII In terms of the variable x = k/t, the size of this 
region ~ i -1 / 2 decays with time. Thus, a diffusive 
boundary layer separates the lower class and the middle 
class. 



equation for the underlying fitness distribution, we ob- 
tain three types of classes: a static lower class, a mobile 
but disperse middle class, and a mobile but "condensed" 
egalitarian class. The population as a whole organizes 
into four types of societies, three of which consist of one 
of these classes, and a hierarchical society in which the 
lower class and the middle class coexist. Two parame- 
ters, the rate of decline and the advancement probability, 
quantify the competition between advancement and de- 
cline. The overall social organization is determined solely 
by these two parameters. 

The fitness distribution has a very different charac- 
ter in each of the classes. In the lower class, this fit- 
ness distribution approaches a steady state. In the mid- 
dle class, the distribution is self-similar in time and cor- 
respondingly the characteristic fitness increases linearly 
with time. Although agents are upwardly mobile, the 
disparities between agents in the middle class also grows 
indefinitely. In the egalitarian class, the fitness distribu- 
tion follows a traveling wave, so that all agents constantly 
advance, but fitness differences between agents remains 
small. 

Much of the richness of the phenomenology is due to 
the fact that the mechanisms for advancement and de- 
cline are fundamentally different. One requires interac- 
tion between agents, while the other is a single-agent pro- 
cess. This dichotomy is reflected by the master equation 
where the decline terms are linear but the advancement 
terms are nonlinear. As a result, there is no detailed bal- 
ance and the dynamics are non-equilibrium in character. 

It should be interesting to use the advancement-decline 
model to analyze real-world data. One natural applica- 
tion is to wealth and income distributions of individuals, 
where both power-law and exponential behavior has been 
observed 0, A related issue is the wealth of na- 

tions. It is well documented that the wealth distribution 
of countries is extremely inequitable, with 60% of the 
world's population producing just 5.6% of the planet's 
gross domestic product (GDP), another 20% producing 
11.7%, and the remaining 20% of the population produc- 
ing 82.7% of the GDP [2!j. The existence of such a large 
underclass corresponds to a large decline rate in our di- 
versity model and it may be worthwhile to understand 
the social mechanisms for such a large decline. Another 
possibility is sports statistics where the winning percent- 
age distribution of teams plays the role of the fitness dis- 
tribution 301. 



VI. CONCLUSIONS 
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